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The magnitudes of the CKM matrix elements Vtd , Vts , and 
Vtb are determined for the first time without any assumptions 
of unitarity. The implications for the unitarity of the CKM 
matrix as a whole are discussed. 
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I. INTRODUCTION 

The relationship between weak and mass eigenstates of 
quarks, assuming there are three generations, is described 
by the Cabibbo-Kobayashi-Maskawa (CKM) matrix [n] 



Vud Vus Vub 

y = I v,^ V,, Kj 

Vtd Vts Vtb 



(1) 



Unless the CKM matrix is assumed to be unitary the bot- 
tom row elements, Vtd, Vts, and Vtb, are undetermined. 
There are, however, a number of measured quantities 
which are sensitive to combinations of these elements. 

We use our recent determination of |Vff,| from elec- 
troweak corrections to Z decays B to "unlock" the bot- 
tom row of the CKM matrix. From a combined fit to data 
from the LEP, SLC, CESR, Tevatron, and other experi- 
ments we determine all the CKM elements, independent 
of unitarity assumptions, and proceed to test unitarity of 
the full CKM matrix for the first time. 



II. CONSTRAINTS ON THE CKM ELEMENTS 

We use experimental results and the corresponding 
theoretical predictions to construct a chisquare which is 
then minimised to obtain all the CKM matrix elements. 
The following sections describe individual chisquare con- 
straints which are ultimately summed to form the global 
chisquare. 



A. Constraint on |Vi6| 

Li a previous paper |Q we described a new method for 
the determination of |Vf6| from electroweak loop correc- 
tions, in particular to the process Z — > hh. We define a 
chisquare Xjj,(mz, ttij, m/f , as, a, Vts) for the agreement 



of the theory with nineteen experimental measurements. 
Full details of the input parameter values, errors, and 
correlations are given in Ref. [EJ . 

B. Constraint on |V(jVtdj 

The product jVj'^Vtdl is constrained by measurements 
of B^ — B^ mixing. Therefore, using these measurements 
and our constraint on |Vt6| enables us to extract \Vtd\- 
The mixing is expressed in terms of the mass difference 
Amd of the CP eigenstates 



Arrid 



r<2 

^^ FlT7'*T7' l2 2 

^iVtbVtdl m^mB^ 



'Jsa 



VbSo (xt) , (2) 



where ^/Bs^fBd = (201 ± 42) MeV ^ comes from a 
combined analysis of lattice calculations, which are com- 
patible with predictions of QCD sum rules, and the 
QCD factor is r]B = 0.55 ± 0.01 ||. The function 
So{xt = ml/m^) is given by H 



S^{.xt) = -^ 



i-Qxt 



(Sx^hixt 



{xt - ly {xt - If 



(3) 



where mt is the running top mass [g| as required for con- 
sistency with rjB- It is related to the pole mass mt as 
measured by the Tevatron, according to 



mt = mt 



_ A as{mt) 

3 TT 



mt 



SGeV. 



(4) 



Given the measurement of Am^d = 0.472 ± 0.018 ps ^ 
we define the chisquare constraint on ll^f^T^dl 



\td^tb 



{Amd{mt,m 



w) 



0.472)' 



{0Al8Amd{mt,mw)) + (0.018)^ 



(5) 



where the theoretical prediction in the numerator de- 
pends on m,t and myy- We vary m,t explicitly in the 
fit; mn/ is varied implicitly since it is a function of the 
other electroweak parameters which vary explicitly [EJ. 
The two terms in the denominator correspond respec- 
tively to the theoretical and experimental uncertainties, 
the former being dominated by the 20% relative error on 

\/BBjBa- 



C. Constraint on \Vt,Vtb\l\Vcb\ 

The rate for the process b —> sj, which involves a loop 
with a virtual top quark, is proportional to | Vf* Vtb P ■ Our 
constraint on |Vff,| therefore enables us to determine \Vts\- 
The theoretical prediction in the Standard Model has 
been calculated to be ifTll 



BR(6 -> S7) 



th. 



(3.28 ±0.33) X 10"^ 



0.976 



(6) 



The element Vcb enters through the normalisation to the 
measurement of BR(B -^ X^eVe) = (10.4 ± 0.4)% |] 
which reduces the otherwise large error due to the b- 
quark mass uncertainty. CLEO has measured BR(6 -^ 
S7) = (2.32 ± 0.57 ± 0.35) x lO"* @, and ALEPH 
has presented a preliminary measurement of BR(6 -^ 
S7) = (3.29 ± 0.71 ± 0.68) x 10-^ @. We aver- 
age these two results, assuming a common uncertainty 
of 0.3 X 10"'* for the theoretical modelling, to obtain 
BR(6 -^ S7)ox. = (2.60 ± 0.59) x 10""*. The chisquare 
constraint on |Vf* Vf6|/| Vcb | is then 



Xts*tb/cb 



(BR(6 ^ S7)th. - BR(6 ^ S7)cx.)' 



^th. 



+ ^i 



(7) 



where cth. and (Tex. are the theoretical and experimental 
errors on BR(& — > 57) respectively. We neglect the ex- 
plicit dependence on rrit and as but we do include their 
3% contribution to the theoretical uncertainty of 10%. 



D. Constraints on other CKM elements 

The most precise constraints on the element \Vud\ are 
from nuclear beta decay. We use the recent measurement 
of \Vud\ from the Chalk River Laboratories [|ll[ to define 
the chisquare constraining \Vud\ 



And 



^(iKdl- 0.9740)7(0.0005)2 



(8) 



We use combined results from measurements of K^^ and 
hyperon decays B to construct the \Vus\ chisquare 



A^ 



IK, 



0.2205)V(0.0018)2 



(9) 



We use the results of experiments measuring the produc- 
tion of charm from (anti-)neutrino interactions with d 
valence quarks Isl to construct the |Vcd| chisquare 



Xcd 



(iKdl -0.224)7(0.016)2. 



(10) 



Measurements of the Des decay D — > Ke^v M are used 
to construct the | Vcs \ chisquare 



where the error is dominated by uncertainties in the 
hadronic form factors. 

The existence oi b —> u transitions, which depend on 
Vub I , is well established from measurements of the end- 
point of the charged lepton energy spectrum in 6 — > 
Xl^v decays. Such analyses yield the ratio |V„hl/|Vcb| ||] 
from which we derive the chisquare 



Auh/cb 



(iKbl/lKftl -0.08)7(0.02)^ 



(12) 



Recently CLEO reported the first measurements of 
BR (B" -> TT-i+v) = (1.8 ± 0.4 ± 0.3 ± 0.2) x lO"* and 
BR (B" -^ p-(.+v) = (2.5 ± 0.4l°;^ ± 0.5) x 10"'' where 
the errors correspond respectively to the statistical, sys- 
tematic, and model-dependent uncertainties |l2|. These 
results yield the chisquare 



Xub 



(iKbl- 0.0033)7(0.0008)2 



(13) 



The element |Vcb| is determined both from the average 
b-lifetime and from B -^ D*lv decays in the limit of zero 
recoil. The resulting chisquare for |Vcf,| is H 



xlb 



(iKbl- 0.040)2/(0.003)2 



(14) 



III. UNITARITY-FREE FIT 

We fit for all the CKM matrix elements without as- 
suming unitarity of the CKM matrix by minimising the 
chisquare, defined as 

x' = xldiVud) + xlsiVus) + xldiVcd) + xlsiVcs) 
+xlb/cbiVub, Vcb) + x7(Kfc) + xlbiVcb) 
+x'td*tbiVtb,Vtd,mt,mw) 

+Xts*tb/cb(.Vts,Vth, Vcb) 



+xibi'm'Z,mt,mH,as, a, Vtb), 



(15) 



where for completeness we have included all the depen- 
dences of the individual chisquare definitions. The 14 
parameters which vary in the fit are mz, rrit, \ogiQ{mH), 
ctsi'nT'z), ct, and the magnitudes of the nine CKM ele- 
ments. The results of the fit are shown in the second 
column of table |l where the errors include the effects 
of experimental and theoretical uncertainties and all cor- 
relations. The chisquare per degree of freedom for the 
unitarity-free fit is 15.1/(28 - 14) = 15.1/14. The proba- 
bility to obtain a value greater than this is 37% indicating 
that the data are in good agreement with the underlying 
model. 

The bottom row CKM elements, \Vtd\, \Vts\, and \Vtb\ 
are determined for the first time, independent of any uni- 
tarity assumptions. 



x2, = (IK. I -1.01)2/(0.18)2 



(11) 



IV. UNITARITY-CONSTRAINED FIT 

We also fit the data imposing CKM unitarity as a con- 
straint. The four parameters which define a three dimen- 
sional unitary matrix can be uniquely determined from 
the moduli of the elements. We choose the parametrisa- 
tion of the CKM matrix advocated by the Particle Data 
Group m which enforces exact unitarity 



Pi 



JSi3 



Vud = C12C13 
Vus = S12C13 

Vub = siae-'^^^ 

Vcd = -S12C23 - Ci2S23Sl3e"' 

Vcs = C12C23 - Si2S23Sl3e**'^ 

Vcb = S23C13 

Vtd = S12S23 - Ci2C23Sl3e*''i^ 

Vts = -C12S23 - Si2C23Sl3e*'''=' 

Vtb = C23C13 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 



where Cy — cos 



lij and Sij = sin ( 



ij for three angles, 6*12, 
013, and 023, which lie in the first quadrant. The phase 
parameter 613 lies in the range < 613 < 27r. The 9 
parameters which vary in the unitarity-constrained fit are 
mz, mt, logio("T-ff), asimz), a, S12, S13, S23, and cos (^13. 
The results of the fit are shown in table | together with 
the corresponding values of the CKM elements. 

The changes in rrit, ttih, and as, compared to the 
unitarity-free fit, are due to their correlation with |Vtf,|, 
as discussed in Ref. [g|- Since only terms of the form 
cos Si3 appear in this analysis, there is a two- fold ambi- 
guity between the regions < (5i3 < tt, as favoured by CP 
non-conservation in the kaon system, and n < 613 < 2tt. 
Unfortunately, the fitted value of cos(5i3 — — 0.261qy4 
has rather large errors; in fact, the lower error corre- 
sponds to the physical limit of cos(5i3 = — 1. 

The chisquare per degree of freedom for the unitarity- 
constrained fit is 25.5/(28 - 9) = 25.5/19. The probabil- 
ity to obtain a value larger than this is 14% which indi- 
cates reasonable consistency of the data and the model, 
albeit with less probability than the unitarity-free fit. 



V. TESTS OF CKM UNITARITY 

Unitarity implies VV"^ = V'^Y = [/, where U is the 
unit matrix, and hence the normality conditions 



Pi 



1^.1 1 



\y^2 



114.'' 



w^ 



2j 



\y^^ 



3jl 



and 



(25) 
(26) 



for 1 = 1,2,3 rows and j — 1,2,3 columns. From only the 
magnitudes of the nine CKM elements it is not possible 
to test orthogonality of pairs of different rows or pairs of 
different columns. Allowing for correlations, the results 
of the six normality tests are 



92 = \Vc 

Pz = \yt 



cd 



td\ 



Kl = \Vud\ 

IK., I 



K2 

K3 : 



\Vudr + \Vus\ 

\Vcs 
\Vts 

\Vcd 
\Vcs 
1^6 



\Vub 



Wu, 



0.997 ±0.001 



\Vcb\ =1.105 ±0.367 
\Vtb\^ =0.623 ±0.315 
\Vtd\'^ =0.999 ±0.007 
WfJ'^ = 1.104 ±0.367 



Wt 



tb\ 



0.629 ±0.316. 



(27) 
(28) 
(29) 
(30) 

(31) 
(32) 



The respective similarities between the three row con- 
straints and the three column constraints reflects the rel- 
atively small values of the off-diagonal elements. Only 
pi is significantly different from unity, being 2.1 standard 
deviations low. Given that there are six, albeit somewhat 
correlated, conditions we do not attach great significance 
to this although its future evolution is of interest. 

The trend is for the p, and Kj to be somewhat lower 
than the expectations of a unitary matrix, as would be 
the case if there were more than three generations. To 
test this possibility with a single number we define the 
quantity 



^CKM 



\/ {PiP2Pz){kiK2H3) 



(33) 



which should be unity for a unitary matrix. We find 
AcKM = 0.69 ± 0.43 where all correlations have been 
taken into account. This is lower than unity but only 
by 0.7 standard deviations, indicating that there is no 
compelling evidence for a deviation from unitarity. 

These results constitute the first tests of unitarity of 
the complete CKM matrix. In particular, p3, ki, k,2, K3, 
and AcKM are measured for the first time. 



VI. OUTLOOK 

The error on \Vtd\ is limited by the large theoreti- 
cal uncertainty on [Bs^f^ ), which is improving only 
slowly with time. In future \Vtd\ will also be deter- 
mined from the process K'^ -^ -k^vD, which is sensi- 
tive to the product |V(^Vts|. Recently, the E787 exper- 
iment observed one candidate event with an estimated 
background of (0.08 ± 0.03) events, from which they de- 
termined BR(/s:+ -^ TT+v v) = (4.2+1;^) x 10"!° @. 
From this they derive, assuming unitarity, the constraint 
0.006 < I Vtd I < 0.060, which is in agreement with our 
more precise unitarity- independent result. Ultimately, 
the measurement of BR(ii''^ — > tt+i/ D) will yield precise 
and theoretically well understood constraints on | V^^ Vts |. 

The uncertainty on |Vts| from the process 6 ^ S7 is 
statistics limited, so some improvement can be expected 
from CLEO and the B factory experiments. The param- 
eter Arris describing B^ — Bs mixing is also sensitive to 
I Vts I . The latest experimental limit is Arus > 10.2 ps~^ 
at the 95% C.L. |^. If Atos were to be measured, its 
interpretation in terms of |Vts| would suffer from the 



large errors on (Ss^/^ ). The ratio Amd/^nis is, how- 
ever, considerably less sensitive to theoretical uncertain- 
ties and its measurement would yield an important con- 
straint on the ratio |T4d/VfsP- Similarly, a future mea- 
surement of the ratio BR{B -^ {p/uj)-f)/BR{B -^ K*-f) 
would constrain \Vtd/Vts\'^. 

Given that LEP has finished running on the Z, the un- 
certainty on I Vtb I from our method of fitting electroweak 
data is unlikely to go below approximately 20% Q . The 
single top quark production rate at hadron colliders is, 
however, also sensitive to \Vtb\- The estimated uncer- 
tainty at the end of the Tevatron Run II is 5|Vtb|/|Vtf,| = 
12% — 19%, where the range reflects the uncertainty of 
the gluon structure functions [|l4| . 

VII. SUMMARY 
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The elements \Vtd\, \Vts\, and \Vtb\, and hence the full 
CKM matrix, are determined for the first time indepen- 
dent of any unitarity assumptions: 



|T/fd|= 0.0113 

\Vts 



-0.0060. 
-0.0029' 



045+°°22- 



1^*61=0.77- 



+0.18 
-0.24- 



(34) 
(35) 
(36) 



We test unitarity and determine that all six CKM nor- 
mality conditions are consistent with unitarity. Four of 
these tests are performed for the first time. 
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TABLE I. Results of the unitarity-free and unitar- 
ity-constrained fits described in the text. Values marked by 
a dagger (f) are not explicitly free in the fit but are derived 
from the other parameters in the same column. 





Free fit 


Gonstrained fit 


mz (GeV) 


91.187 ±0.002 


91.187 ±0.002 


mt (GeV) 


174.2 ± 5.7 


172.4 ± 5.3 


logio("^H/GeV) 


9 -|r + 0.30 


2.03l°;f, 


Qs(mz) 


0.1171 ±0.0025 


0.1188 ±0.0021 


a~'^{mz) 


128.913 ±0.092 


128.904 ± 0.092 


\Vud\ 


0.9740 ± 0.0005 


0.9748 ± 0.0003 f 


|v;,| 


0.2205 ±0.0018 


0.2230 ± 0.0014 t 


jKil 


0.00325 ± 0.00058 


0.0032 ± 0.0006 t 


\Vcd\ 


0.224 ±0.016 


0.2228 ± 0.0014 f 


\V,s\ 


1.01 ±0.18 


0.9741 ± 0.0003 t 


\Vcb\ 


0.0401 ± 0.0029 


0.0397 ± 0.0029 f 


\Vt,\ 


omistlZll 


0.0093 ± 0.0018 t 


\Vts\ 


0.045lHfo 


0.0387 ± 0.0029 f 


\Vti\ 


0.77tlM 


0.9992 ± 0.0001 t 


Sl2 


— 


0.2230 ± 0.0014 


Sl3 


- 


0.0032 ± 0.0006 


S23 


— 


0.0396 ± 0.0030 


cos 5l3 


- 


-0.26+°;?^ 



